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Abstract 

In a realistic theory with large extra dimensions and TeV scale gravity, 
there are often other particles hving in the extra dimensions in addition to 
gravitons. A well motivated candidate is the flavon in a flavor theory which 
explains the smallness of the fermion masses by the large volume of the extra 
dimensions. The flavon interactions involve the Higgs fleld and can therefore 
give rise to new mechanisms of Higgs boson production. We study these new 
Higgs signals, including both real flavon emission and virtual flavon exchange. 
They could provide powerful probes of the string scale which are comparable 
to the extensively studied graviton production, if the flavons and the graviton 
feel the same number of extra dimensions. If the flavons only live in a sub- 
space of the extra dimensions, they can yield even stronger bounds and may 
generate spectacular displaced vertex signals. Our study also applies to a four 
dimensional flavor theory with light flavons as well as to any theory with bulk 
scalar flelds coupled to the Higgs and fermions. 
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1 Introduction 



The Standard Model (SM) is very successful in explaining most precision data from current 
collider experiments, but still we have not found (or tested) all ingredients of the SM. The 
elusive Higgs particle is yet to be discovered, and its interactions with the SM fermions 
and gauge bosons will need to be measured before we close the final chapter in the quest 
for the correct theory of particle physics at the weak scale. 

Given the success of the Standard Model, there are still several unresolved puzzles. An 
important one is related to flavor, i.e. the size of the fermion masses and the generation 
of the CKM matrix. In the Standard Model, the fermion masses arise from the Yukawa 
interactions between the fermions and the Higgs. We do not understand why all fermion 
Yukawa couplings other than the top Yukawa coupling are so small or why the CKM 
matrix is approximately diagonal. The smallness of the fermion Yukawa couplings also 
leads to the suppression of direct s-channel Higgs production, making the Higgs discovery 
a challenging task. 

A natural explanation for the smallness of the fermion Yukawa couplings is that there is 
an enlarged (flavor) symmetry under which the direct couplings between the light fermions 
and the Higgs are forbidden in the symmetry limit. Only after the flavor symmetry is 
broken by the vacuum expectation value (vev) of some scalar fiavon field {(f), the Yukawa 
couplings of light fermions are generated by some higher dimensional operators which are 
invariant under the flavor symmetry, such as 

§;^Hhh, (1.1) 

with ?/j ~ 1. The Yukawa coupling is given hj yj {ip) /Mp and is small if (y?) -C Mp. 
For a more general flavor theory, there can be several yj's and the Yukawa couplings can 
be suppressed by higher powers of {if) jMp. The higher dimensional operators like (|1.1|) 
can arise from gravitational scale physics, in which case Mp will be the Planck scale, or 
from integrating out some vector-like Froggatt-Nielsen fields [|1|, in which case Mp will be 
the mass of the Froggatt-Nielsen fields. The scale of the flavor physics is not determined 
by ( |1 . 1| ) as long as the ratio ((/?) jMp produces the necessary fermion mass hierarchy. 
However, the flavor symmetry breaking scale in the usual four dimensional theory is 
typically quite high. If the flavor symmetry is a continuous global symmetry, there will 
be one or more massless Goldstone bosons, in which case there are strong constraints from 
experiments and astrophysics requiring the symmetry breaking scale to be greater than 
~ 10^ GeV. In the cases of a gauge symmetry or a discrete symmetry, the gauge bosons 
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and the flavons will induce flavor- changing effects, depending on the symmetry and its 
breaking pattern. On the other hand, the absence of small couplings in the fundamental 
theory implies that the masses of the fiavons are on the order of (or smaller than) their 
vevs. Since the masses are constrained from experiment (the non-observation of new 



the weak scale, the fiavon vevs themselves must lie in that range too. Then, to get a small 
ratio of (y?) /Mp, one has to consider values of Mp much larger than the weak scale. 

Another unsettling issue with the SM is the hierarchy problem. A few proposals 
(supersymmetry, technicolor) have been around for a long time, but the most radical 
approach 0] was put forth only recently. The idea is to introduce extra large compact 
spatial dimensions to bring the fundamental physics scale down to a few TeV, thus ame- 
liorating the hierarchy problem. Without the large hierarchy between the fundamental 
Planck scale and the weak scale, any flavor physics must occur close to the weak scale. In 
fact, without the grand "desert" to protect the SM from various problems, such as proton 
decay and large flavor-changing effects, it is desirable to have a large (discrete) flavor sym- 
metry to prevent these disastrous large effects 0. The hierarchy of the fermion masses 
in this framework can originate from the large size of the extra dimensions instead. For 
example, the flavor symmetry may be maximally broken on one or several distant branes 
in the extra dimensions. The breaking is then transmitted to the brane where the SM 
fields reside by some bulk fields with the same quantum numbers as the flavor symmetry 
breaking fields on the distant branes. Depending on the distances and the masses of the 
bulk fields, the symmetry breaking effects on our brane may be suppressed by powers of or 
exponentially with the distance, which explains the smallness of the Yukawa couplings ^ . 
Specifically, the interactions on our brane look like 



where x is the bulk messenger field (with mass dimension 1 +5/2, 6 is the number of the 
extra dimensions accessible to x)i U represents the coordinates of the extra dimensions, 
and Ms is the fundamental scale of the higher- dimensional theory, possibly the string scale. 
X has a nonzero {y dependent) vev induced from the distant brane. After propagating 



framework, the smallness of {x{x,y = 0)) comes from the distance suppression and does 
not require the mass of x to be of order its vev (x), in contrast to a four dimensional 
theory. 



gauge bosons and compliance with rare flavor- changing processes 0, ^) to be well above 




(1.2) 



the distance between the two branes, on our brane x has a vev {x{y = 0)) <^ M]+'5/2 and 
as a result, the corresponding Yukawa coupling Xj = y'j {x{y = 0)) /M}^^^'^ <^ 1. In this 
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The string scale is expected to be close to the Planck scale in the higher dimensional 
theory. For n extra dimensions accessible to the gravitons (n > 5) , the 4 + n dimensional 
Planck scale M^, and reduced Planck scale M^, are related to the 4 dimensional Newton's 
constant Gat by 

M2+"V; = M2+"K/(2vr)" = {A-kGmV = ^1, (1-3) 

where Vn is the n dimensional volume of the extra dimensions {Vn = if all extra 
dimensions have the same length L). To avoid retaining too many similar scales and to 
simplify the notation, we will use the 4 + n dimensional reduced Planck scale for most of 
the paper and define 

M = M,. (1.4) 
The interaction ( [1.2|) can be rewritten in terms of M, 

J d'xhfnHxix, y = 0). (1.5) 

Remember that gravitons have order one couplings in units of the reduced Planck scale. 
We will assume ~ 1 for our study. 

The bulk messenger field x (which from now on we shall call a "fiavon" , in analogy to 
the case of a four dimensional flavor theory) can live in a subspace of 4 + 5 dimensions 
with 6 < n. From a 4 dimensional point of view, the x fi^ld is Fourier decomposed into a 
Kaluza-Klein (KK) tower of states, 

x(l/ = 0) = ^EXn, (1.6) 



where Vs is the volume of 5 extra dimensions (V^ = for equal length), and n's are 
integer vectors in the 6 dimensions. The operator (|1.5|) can then be written as 

EhhHxn, (1.7) 



where we define 



Mf= M^+^/V//', (1.8) 



so that the interaction has the same form as ( |1 . 1|) . If x propagates in the same extra 
dimensions as gravity, i.e. if 6 = n, Mp can be identified with the 4-dimensional reduced 
Planck scale Mp = (47rG'Af)~^/^. However, Mp can in principle lie anywhere between M 
and Mp. For example, if all extra dimensions have the same length L, we find Mp = 

MiMp/Mf/'' = Mp{M/MpY-^/''. 



3 



Although each interaction ( |1.7| ) is highly suppressed by the large scale Mp, the large 
number of KK states compensates for the suppression. In the end of the day, the cross- 
sections will only depend of M and 5 and not on Mp or L, which we only introduce for 
notational clarity. 

For very light bulk messengers, i.e. "fiavons", (with zero mode mass mo < 1 GeV), 
there are various constraints from astrophysics and flavor violating effects pushing the 
lower bound of M to be far beyond the weak scale for 5 < 4 All these constraints 
disappear if mo is greater than a few GeV, because the x KK modes would be too heavy 
to be produced in these processes. In that case, M is only constrained by the graviton 
KK state interactions and can be ~ TeV for n > 2. For such a low M, the interaction 
( |1.5|) can be important at high energy colliders and give rise to new Higgs signals in future 
experiments. 

The possibility of x living in fewer dimensions than gravity may result in some advan- 
tages of collider probes of the x fields compared to the more model independent signals 
from KK gravitons. First, the resulting limits on the gravitational scale in collider exper- 
iments are stronger for smaller number of extra dimensions 0. However, the existing 
lower bounds on M are also stronger for smaller n = 1 is ruled out, M^, > 50 TeV from 
SN1987A and M^, > 110 TeV from cosmic diffuse gamma ray radiation for n = 2 [^). 
As mentioned above, these constraints do not apply to x if X is heavier than a few GeV. 
It is possible that n is large, yet 6 is small, so that M is quite low and we have larger 
collider signals for x rather than the KK states of the graviton. Second, because Mp can 
be smaller than Mp for 6 < n, there is a range of Mp where the KK modes of x fields 
can have easily detectable macroscopic decay lengths inside the detector. 

The interactions between different fermions and x's depend on the flavor symmetry 
and how it is broken. A large flavor symmetry is favored to avoid the flavor violation 
problem coming from higher dimensional operators generated at the string scale Mg. A 
straightforward example for the flavor symmetry is (a large discrete subgroup of) U(3)^ ||^, 



10|. The X fields transform as (3, 3) under the product of two U(3)'s. The vevs of x's 



break U(3)^ and generate the SM fermion Yukawa couplings. It is known that the flavor 
changing effects are safe if the SM fermion Yukawa couplings are the only source of U(3)^ 



breaking [11]. Additional flavor violating effects may be induced from the x interaction 
away from our brane but are suppressed by distances . It provides a working example 
and for the most part of this paper we shall work under this assumption, that is, we will 
assume that there is one and only one x field for each entry of the SM fermion Yukawa 
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Table 1: Current limits on the scale M (in GeV) from real graviton emission 
processes, for different number n of extra dimensions. 



matrices, e.g., HqiqjXqiqj HliCjXii^y For a more general flavor theory, there may be x 
fields which couple to more than one combinations of fermions, and the Yukawa couplings 
may come from interactions containing higher powers of x fields or products of different 

To summarize, the scenario that we are considering here is completely determined by 
just a few independent parameters (assuming that all extra dimensions have the same 
size): 

{M, mQ,mh,5,n} . 

Mp and L are derived quantities: Mp = M{Mp/MY/'' and L = M-^{Mp/Mf/''. The 
allowed range of M and n is constrained by the astrophysical and collider bounds on 
gravity in large extra dimensions. In Table |I] we list the current limits on the scale M from 
real graviton emission at LEP and the Tevatron, translated from the studies of several 
groups P, 1^, [1^, [13|. We do not include bounds from processes with virtual gravition 
exchange, where the appearance of an explicit ultraviolet cutoff in the calculation makes 
it difficult to interpret the limit on M. As we mentioned earlier, tuq has to be larger than 
a few GeV, while m/j has to satisfy the experimental limit from the Higgs search. 

The plan of the paper is as follows. In Section we first check whether the presence 
of the additional operator ( |1.7|) can affect the conventional Higgs searches at LEP and 
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the Tevatron. In the next two sections we explore the new discovery signatures offered by 
the couphng ( |L7D . In Section || we concentrate on the direct production of KK states Xn, 
leading to rather unorthodox Higgs signals associated with missing energy or displaced 
vertices. Then in Section § we consider new Higgs production processes with a virtual 
exchange of KK states Xn- In all cases, we estimate the discovery potential of both LEP 
and the Tevatron in the new channels, and also comment on the possible reach of LHC 
or NLC. We summarize and draw our conclusions in Section ^. 



2 Conventional Higgs searches 



The traditional searches [|T^, [T5[ look for the SM Higgs in associated production // — »• 
Wh,Zh and gluon fusion gg h. A two-Higgs doublet model, e.g., the MSSM, also 
presents the opportunity to look for pp — > hbb at large values of tan/?. The specific 
final state topology depends on the Higgs mass m/j. For rrih ^ 140 GeV, the Higgs decays 
predominantly to bb pairs, while for nih > 140 GeV, the decays to W^W~ pairs dominate. 
All of these searches will be significantly affected if the higher dimensional couplings to 
fiavons can change dramatically the branching ratios of the Higgs to bb or W^W~ . Indeed, 
r(/i — > bb) is quite small, which has already prompted suggestions to look for other (e.g. 



invisible |]T6|, [T^) Higgs decay modes for low Higgs masses. At the same time, the Higgs 
search is the cornerstone of the Tevatron Run II physics program, peaking in sensitivity 
for a range of Higgs masses just beyond the LEP-II reach. As for somewhat heavier 
Higgs masses, nih ^ 140 GeV, the width V{h W^W~) is already rather large, so the 
branching ratios to gauge boson pairs are unlikely to be affected by any new physics. 

In order to estimate the effect of the operator ( |1.7|) on the Higgs branching ratios, we 
have to compare the width V{h bb) to T{h xJj) or T{h — ^ x^'^^')- To leading order 
in the fermion mass m/, one has 

where A^^^ is the number of colors of / and M^y is the W mass. In the presence of ( p,.5| ), 
and for mo < rrih, three-body decays to a fermion pair and the corresponding fiavon are 
also allowed. For a single KK mode of x? 

— A^ vt^Tfx^ r 1 'I 
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where /i = m\/m\ and rrin is the Xn mass. Summing over the KK modes we have 

^{h^Xff) = dmlp{m^)T{h^Xnff) (2.11) 

J mi 



1 N,y} mi+^ 

(47r)V2r(5/2) 2567r3 M^+^ 



where p{m-n) is the density function of the KK modes given in the Appendix, and we have 
defined /iq = ml /ml and 

Gsipo) = £ dfiifi - p,f~^y^(^^(l + 9/x - 9/x2 - + + iog^|_ (2.12) 

In the hmit of mo 0, the G5(0) integral can be performed analytically and the result is 

128 

= 5(2 + 5)^(4 + 5)2(6 + 5)' ^^-^^^ 

which ranges from 8.1 x 10"^ for 5 = 1 to 2.8 x 10~^ for 5 = 6. The ratio of the widths 
of the new h xff channel to h ^ bb is then given by 

.2 / »^ \ 2 / \ 5+2 



T{h bb) 3 87r2 g"^ \ m^ ) \M ) (47r)'5/2r(5/2) 

" (47r)V2r(5/2) (?2 l^Myl 



:2.14) 



We show this ratio in Fig. |l] for fixed m^ = 120 GeV, mg = 10 GeV, = 3 and yf = l 

(for larger mo, it will be even smaller). In the figure we also mark the typical limits 

from Table |l| on the scale M for various values of n. We see that the branching ratio to 

flavons is typically orders of magnitude smaller than the branching to bb, and is therefore 

unlikely to have any impact on the traditional Higgs searches^ The only exception is 

the case of 5 = 1, where it can be as large as a few percent. One should not forget that 

the width being compared to r(/i ^ bb) in Fig. |l| is for a single flavon decay mode. In 

principle, one should sum over all fiavon flavors, which give similar final state signatures, 

e.g. Higgs decays to two jets plus missing energy can be obtained from decays to any 

Xqiqj-i = u,d,s,c. This will enhance the branching fraction, but one loses the b- 

tagging option, which is crucial for reducing the background in a hadronic environment. 

On the other hand, there is also the possibility for decays to leptonsQ via h ii'xee', 

^Our conclusion in this respect differs from tliat of Rcf. [||, where the numerical factor 
Gs{^io)/ [87r2(47r)*/2r((5/2)] was not included. 

^From this point on we use the term "lepton" {£) to denote a muon or electron only. 



7 




200 300 400 500 



M (GeV) 

Figure 1: The ratio of the Higgs widths r{h — > xff)/^{h bb), as a function of 

M, for different values of S, and for fixed rrih = 120 GcV, mo = 10 GeV, Nc = 3 
and Df = 1. The arrows are to remind the reader of the typical hmits on M from 
graviton emission processes. 

albeit with a much smaller branching fraction (by a factor of 7.3) than h — > jjXqq- The 
final state is indistinguishable from the SM decay h ^ tt ^ ii/^i'Pi'UrpT, which already 
has a branching ratio of about 0.6%, much larger than what we typically would expect for 
h — > ii'xw- Since the tau decays are uncorrelated, the different fiavor decays {£ ^ I') are 
not particularly distinctive either, besides, the major backgrounds from tt and W^W~ 
half of the time also yield different flavor leptons. 

In conclusion of this section, we see no easy way to discover or conflrm this scenario 
by looking at the decay modes of the Higgs particle. One (rather exotic) option for that 
would be to build a muon collider at the Higgs mass and then measure precisely the tau 
decay modes of the Higgs, looking for nonuniversahties. Even then, on has to hope that 
the flavon couplings to leptons exhibit some kind of nonuniversality, otherwise the Higgs 
to tau decays would look completely normal, with a somewhat higher branching fraction, 
which in turn may be attributed to a number of other sources of new physics. At a muon 
collider, it will be better to look for the h — > jjx modes, which are more abundant than 
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Table 2: The decay length cr in meters of a KK flavon of mass 1 GeV for 
M = 1 TeV and various number of extra dimensions n and 6. 

h iix anyway. Again, they will have to be distinguished from h ^ tt ^ H^t^t) which 
certainly seems possible. First, the missing energy in the decays of a highly boosted tau 
is correlated with the jet direction, and second, hard tau jets are quite different from light 
quark jets jTSl . 



3 Direct Flavon production 



In this section we shall consider the production of a real flavon in association with a real 
Higgs in both e^e~ and pp collisions. The resulting signatures will depend on the fate of 
the flavon KK states, hence we discuss their lifetime first. 



r(Xn-//) = T^^TTT-n. (3.15) 



The width of a real Xn decaying back to the corresponding fermion pair is given by 

Nc yj 

Depending on Mp (or equivalently, the size of the extra dimensional space in which x 
can travel), Xn can decay promptly at the interaction vertex, have macroscopic decay 
length inside the detector, or decay well outside the detector. If one assumes that all 
extra dimensions have the same size L, then the decay length of the KK flavon state Xn 
can be expressed in terms of 5 and n, 

e.(.,^(^)='--'(i^):o--., 

and is shown in Table ^ for various n and 5. 

For prompt decay, the signal is the same as the virtual x exchange which will be 
discussed in the next section. If x-n has a macroscopic decay length inside the detector. 
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we will have the spectacular signal of displaced vertices . Displaced vertices can appear 
in several versions of the supersymmetric extensions of the Standard Model. An example 
which gives similar signatures is the low supersymmetry breaking scale model with neutral 
wino as the NLSP For a range of values for the supersymmetry breaking scale, the 
NLSP has a delayed decay to a Z and a gravitino, with Z decaying subsequently to a pair 
of leptons or jets. Note that the decay length also depends on the mass of the KK mode 
and we have a continuous distribution of KK mode masses, so the range of Mp which can 
give rise to the displaced vertex signature is not too narrow.^ 

Finally, if the lifetime of Xn is large enough, it will escape the detector. This will give 
rise to an interesting ^ signature with real Higgs bosons, somewhat similar to Zh — > vvh. 
The difference now is that the missing mass is not characteristic of the Z mass, but rather 
displays a continuous spectrum. If we replace the Higgs by its vev in the interaction and 
attach a photon or a gluon, we can also get signals similar to KK graviton production 
such as 7+ ^ or j + Ip. Since they have already been studied for the case of graviton 
production 0] and do not have Higgs associated with it, in the rest of this section we 
will instead concentrate on the novel h + J/^t signals. 

At both lepton and hadron colliders the production of hxn will lead predominantly to 
a hh^T final state for < 140 GeV, and WW^t or ZZ^t for ruh > 140 GeV. 

The (parton level) cross-section for producing a single fiavon species of mass mn is 
given by 

where 

\{x,y,z) = + + — 2xy — 2xz — 2yz (3.18) 

and is the center of mass energy of the collision. To get the total cross-section (ytot{.s), 
one needs to integrate over the available mass range for the KK states {mo, a/s — rrih}: 

0-totis) = / dmn2mnp{mn)an{mn,s) (3.19) 

J mo 

- 64^^M^ (47r)V2r(|) Lo ^^"2m„(m„ m,) — 

In the remaining of this section we shall discuss the case of lepton and hadron colliders 

separately. In the latter case, the cross-section ( |3.19|) in addition has to be integrated over 

■^In principle the distribution of the decay lengths will be different from that of a particle of a single 
mass. However, with a low statistics sample, it will be very difficult to distinguish these two cases 
experimentally. 
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Figure 2: Signal cross-section for hx production at LEP with y/s = 200 GeV, for 
nih = 120 GeV and different values of 6 and M. 



the respective parton distribution functions f '^{x), e.g. at the Tevatron we have 
with Ecu = 2 TeV now being the total center of mass beam energy. 



1 rl r^x^X2Ecm-mh 

dxi / dx2 / dmn2mnpimn)nxi)fP{x2)animn,XiX2Ec^) 

Jo Jmo 

(3.20) 



3.1 Lepton colliders 

For rrih + mo ^ 200 GeV, KK modes of the flavon can in principle be produced at LEP. 
The total cross-section ( p.l9| ) for e'^e" hx with ^/s = 200 GeV is shown in Fig. |^, for 
ruh = 120 GeV and several values of 6 and M. For a different M, the cross-section can 
be easily obtained from the scaling relation a oc M~^^^^\ The cross-section maxes out at 
small values of mo, when there are more KK states kinematically available. From Fig. ^ 
we see that if 5 = n, the typical maximal values of the cross-section are around 10 fb, 
in which case only one or two events can be expected at LEP. Of course, the bounds on 
M from gravity are weaker if one considers cases of 6 < n, then LEP might be able to 
explore the region of very low values of M and small 6. Given the current limits on the 
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100 200 300 400 500 

(GeV) 

Figure 3: Contour plot of the total cross-section for real flavon emission at the NLC 
with ^/s = 500 GeV, for = 1, 5 = 2 and M = 1 TeV. The shaded area is not 
kinematically accessible at this center-of-mass energy. 

Higgs mass, one obviously has to go to a higher energy machine in order to start to see an 
appreciable signal. We therefore turn to discuss the NLC case in a little bit more detail 
next. 

In Fig. ^ we show a contour plot of the total integrated cross-section at an NLC with 
^/s = 500 GeV. We fix i/f = 1, 6 = 2, M = 1 TeV and vary the Higgs mass mh and 
the mass tuq of the fiavon zero mode. The cross-section shown in Fig. |^ can be easily 
rescaled for different values of ?// or M, recall that it scales like y'j/M'^^^ . With higher 
energy and larger integrated luminosity, the fiavon production cross-section at the NLC 
can provide a much more statistically significant signal sample. With 50 fb~^ of data 
for the parameters in Fig. ^ one will have a few hundred signal events available. We 
will then focus our discussion on the question of cuts optimization and increasing the 
signal-to-background ratio. 

The search strategy at the NLC will depend on the Higgs mass. For heavy ( > 140 
GeV) Higgs masses, the signal will be WW similar to chargino pair production with 
charginos decaying to real 14^'s and the lightest neutralino. In that case, it is best to 
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reconstruct the W^s from the four jet final state [PT|, where the only non-negligible back- 
ground is WW Z with Z — > z/z/, which is on the order of 10 fb after taking into 
account the Z branching ratio into neutrinos. With cleverly designed cuts |^3|, it is easy 
to see an excess over the background for a signal of 0{1 fb), which can then be used to 
put constraints on mo, m/i, and M. 

If, however, the Higgs mass is below 140 GeV, the search becomes more challenging, 
since the backgrounds to the hb^ final state are much larger. The dominant background 
contributions come from bbveVe-, bbv^v^ and hhvT-Vr final states. Including the SM diagrams 
with a Higgs into the background, one has a total of 23, 11 and 11 diagrams, respectively. 
We show the diagrams for the process e^e" ^Vei^e in Fig- i- The diagrams for the 
processes e^e" —>■ hbv^^v^^ and e^e~ — * hhu-rVT- are obtained by obvious replacements and 
dropping the diagrams containing a W boson. Notice that most often the two b-jets are 
not coming from a Higgs, so one can use a {irih dependent) invariant mass cut on the b- 
jets. The optimal mass window will depend on the features of the particular detector, for 
our study here we shall take \mbb — mh\ < 10 GeV. In addition, the availability of highly 
polarized beams at the NLC can substantially suppress the background (which necessarily 
involves the production of W or Z bosons), since W^s only couple to left-handed leptons. 
Finally, the ZZ background with one Z decaying invisibly can be further suppressed by 
requiring that the missing mass il/f in the interaction is not near Mz- The remaining 
WWZ background is small and can be further reduced by a lepton veto. 

We used the COMPHEP event generator to simulate the bbuu backgrounds (with 
separate runs for electron and muon/tau neutrinos) and in Fig. |^ show the resulting 
missing mass ]]/[ distribution at the NLC with ^/s = 500 GeV, where we assume a 90% 
right-polarized electron beam. We also present signal distributions for = 120 GeV, 
mo = 120 GeV and several values of 5 and M: (a) 5 = 1, M = 1940 GeV, (b) 5 = 2, 
M = 858 GeV, (b) 5 = 3, M = 514 GeV and (c) 5 = 4, M = 360 GeV, where we have 
picked M so that the total signal cross-section is always 20 fb. We choose to do this since 
our objective here is to compare the shape of the signal distributions to the shape of the 
background. 

We see from Fig. ^ that the shape of the missing mass distribution for the signal 
depends on the number of extra dimensions S. In particular, for 5 = 1 it exhibits a narrow 
asymmetric peak near mo, for 5 = 2 it is rather flat, while for larger 6 it starts to peak 
closer to the largest kinematically allowed values and therefore becomes very similar to 
the background distribution. (This behavior is easy to understand in terms of the density 
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Figure 4: The diagrams for the process e^e hbUf-Ve- 
bosons apphes for e~''e~ — > bbv^Uf^ and e~''e~ — > bbi^ri^T- 



diagr. 23 

The subset without W 



function - see eq. ( |A.27| ) from the Appendix.) Because of the finite detector resolution, 
the narrow one-sided peak for 5=1 will be smeared and the flavon will look like an 
ordinary heavy, invisibly decaying neutral particle, so if a signal is seen, an interpretation 
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Figure 5: Missing mass distribution of signal and e^e~ hhviVi background at 
the NLC with center-of-mass energy y/s = 500 GeV. The signal distributions are 
shown for nrih = 120 GeV, mo = 120 GeV, and several values of 5 and M: (a) 5 = 1, 
M = 1940 GeV, (b) 5 = 2, M = 858 GeV and (c) 5 = 3, M = 514 GeV, (d) 5 = 4, 
M = 360 GeV. All lines are normalized to the correct total cross-section. We have 
assumed 90% polarization and imposed the cut |m;,fe — m^l < 10 GeV. We show the 
bbveVe background separately from hbv^v^ + bbvri^T (shaded). 



in terms of the scenario considered here will not be straightforward. Notice that there is 
very small background with i^f < Mz, so by far the easiest case to observe will be 5 = 1 
with mo < Mz, since signal and background are very well separated. 



3.2 Hadron colliders 

In Run II of the Tevatron one can look for a bb^T final state, which is one of the standard 



Higgs search channels [25|, but without mass reconstruction may also be indicative of 



several other sources of new physics: bottom squark production ||26|, stop production 
with Higgsino LSP leptoquarks etc. We can adapt the existing Run II Higgs 
analysis of Zh —>■ bb^^ for our purposes, keeping in mind that further optimization is 
possible, e.g. by adjusting the cut and the px cuts on the b-jets. 
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Figure 6: Plot of the signal cross-section for real flavon emission at the Tevatron 
for mh = 120 GeV, ttiq = 100 GeV, M = 500 GeV and several values of 5. 



In order to make an estimate of the signal cross-section at the Tevatron, let us first 
switch the order of integration in (p.20 ): 



crtotiEcu) = / dmn2mnp{mn) dxi , dx2r{xi)f'^{x2)a^{mn,xiX2EcM) 

J mo Jx{mn) J£i£lai 

(3.21) 

where x(mn) = {rrin + / E^j^. Then, as long as we are interested in the total signal 
cross-section only, and not in the resulting kinematic distributions, we can use COM- 
PHEP, for example, to first compute the differential cross-section at a given rrin 

an(mn,EcM) = / = [ dxi j 6X2^X1) f'P{x2)(Tn{m^,XiX2ElM) 

(3.22) 

and then numerically integrate over uin- 

In Fig. P we show a plot of the differential signal cross-section crn(7Tin, 2 TeV) before 
cuts. Here we include the signal contribution only from real fiavon emission and not from 
Zh. Notice that depending on the flavor of the quarks participating in the hard scattering, 
different species of fiavons will be produced. However, we assume that they all escape the 
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detector and give rise to I^t- For the plot we have made the simphfying assumption that 
the couphngs yj are universal and that the mass of the zero mode mo is the same for all 
species of flavons. 

In Fig. ^ we observe the typical dependence of the mass distribution on the number of 
extra dimensions 5, already familiar from Fig. ^ We see that for 5 > 2 the distributions 
peak at rather large values of m-n of around 300-400 GeV. Notice, however, that although 
the actual missing energy in those events will be very large (typically on the order of a few 
hundred GeV), the distribution of the measured will be much softer. The measured 

in the detector is determined by the energy of the visible products, in our case the two 
b-jets, which is not much above m/j/2. Nevertheless, we remind the reader that the Higgs 



search in Run lib will be sensitive to raw bb^x cross-sections on the order of 26 fb |25 
We then see from Fig. |] that flavon searches in the same channel may yield very strong 
bounds on M, for example, we can see from the plot that M < 1000, 500, 300, 200 GeV 
for 5 = 1, 2, 3, 4 respectively can be easily ruled out. 

At a pp collision machine like LHC, the initial state anti-fermion must come from 
a sea quark, which will suppress the production cross-section. At the LHC it will be 
very hard to observe the bb^T final state, instead one typically has to utilize the decays 
fi ^ ryry aud/or k ZZ . Tfic possibility of producing an additional invisible fiavon state 
in the interaction gives rise to unique clean signatures like (for low Higgs masses) 

or llll^T (for higher Higgs masses). 

4 Virtual Flavon exchange 

Throughout this section we shall assume that each fiavon only couples to the corresponding 
pair of fermions. Then, virtual fiavon exchange at an e"'"e~ machine can only give back 
final states with electrons and positrons, while fiavon exchange at a hadron collider will 
yield signatures with jets. 

Virtual fiavons can mediate processes with 2-fermion, 4-fermion or 6-fermion final 
states. The corresponding Feynman diagrams for an e"'"e~ lepton collider are shown in 
Figs. 0, I and ^, respectively (for a hadron collider, just replace e ^ q). 

For 6 > 2, the virtual x exchange can be represented by a local effective dimension 
8 operator ( [A.29| ) as shown in the Appendix. In Table |^ we list the signal cross-sections 



(in fb) from this effective operator at different colliders for A^ = 1 TeV and rrih = 120 
GeV. In the case of a real Higgs production, we only consider decays to bb final states. 
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Figure 7: The diagrams for the process e"^e — > x* ^ 




Final 

state 


LEP 


NLC 


FNAL 


LHC 


200 GeV 


500 GeV 


2 TeV 


14 TeV 


// 


0.95 


5.91 


12.9 


373 


ffh 


9 X 10-5 


6 X 10-2 


2.8 X 10-2 


60 


ffhh 




6.5 X 10-^ 


2.9 X 10-^ 


2.7 



Table 3: Signal cross-sections in fb for the 2-fermion, 4-fermion {ffbb) and 
6-fermion [ffbbbb) final states, at various colliders (with their center of mass 
energy shown) and for = 1 TeV and rrih = 120 GeV. 
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whose branching fraction was computed with the HDECAY program |^. If Ac = M 
and yj = 1, then Ar = 1 TeV corresponds to M = 400,237, 145 GeV for 5 = 3, 4, 5, 
respectively. These cross-sections can be readily rescaled for different values of A7-, recall 
that a ~ A;^^. 

We see from Table ^ that the signal cross- sections are rather small, and just like in the 
case of virtual graviton exchanges, one will have to look for interference effects instead, 
since those are only proportional to A;^^. 

For 6 < 2, the effective interaction still depends on the momentum flow along the x 
propagator. In the case of 5 = 2, there is only a mild logarithmic dependence on the 
momentum and the cutoff. The case of 5 = 1 is more interesting. The integral over 
KK states is convergent, therefore independent of the cutoff. For p'^ > {p is the 
momentum of the x field), the effective interaction is imaginary, corresponding to real x 
production. Whether it gives rise to the same final state particles or not, depends on the 
X lifetime. Real x production does not interfere with other virtual (off-shell) processes. 
The singularity at = itlq will be smoothed out by the decay width of the x fields and 
the detector resolution. Nevertheless, the peak at = provides an extra handle of 
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the 6 = 1 case. 

In the rest of this Section we shall discuss lepton and hadron colliders separately. 
4.1 Lepton colliders 

The virtual scalar exchanges in Fig. |^ will modify the angular distribution of the leptons 



in Bhabha scattering and may be detected at LEP. In Fig. [T^ we show the differential 
cross-section da /d cos 6 (in fb) at LEP for At = 750 GeV, where 6 is the scattering angle 
of the electron with respect to the direction of the electron beam. The shaded histogram 
is just the background contribution, while the solid histogram includes both signal and 
background. There is a visible excess of backscattered electrons, which is typical of the 
(isotropic) scalar exchange. In order to extract the LEP limits on At from Bhabha 



scattering, we can apply the latest bounds on the relevant compositeness scale Alr ||30 
We have to equate the coefficient of the contact term 



-^{eL'j^eL){eR-ff,eR) ~ --r2-(^LeR){eReL) (4.23) 
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100 200 300 400 



mn (GeV) 

Figure 11: Invariant mass distribution of the e^e^ pair in e^e^hP events for both 
signal (solid) and background (shaded) at the NLC with y/s = 500 GeV, and for 
= 600 GeV and mu = 120 GeV. 



to the coefficient of the effective operator ( |A.29| ). The typical bounds on Alr of 4.5 TeV 
for ~ 47r then translate into 

1/2 



At> (^^1 - 525 GeV. (4.24) 



We should point out that in the case of 6 = 1 and mo not near Mz, one may use in 
addition the peak in the invariant mass of the leptons. 

We next discuss the 4-fermion channel, which is similar to the existing Zh e^e~hb 



searches |T^. Here the invariant mass me+e- of the e^e~ pair is not constrained to be 
close to Mzi so one might think that me+e- will be a useful variable to cut on. In Fig. |Tl] 
we show the rrie+e- distribution in e~^e~h^ events for signal plus background (solid) or 
just the background (shaded) at an NLC with y/s = 500 GeV, and for At = 600 GeV and 
rufi = 120 GeV. We see that rrie+e- does not allow good discrimination between signal 
and background, since the two distributions are very similar. Above the Z-peak, the 
background is dominated by Z-fusion, just like W^-fusion was the main background for 
the bb^T channel in Section p.l| . Unlike ly-fusion, here we have the additional advantage 
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Figure 12: The same as Fig 11, but for the angular distribution da/d9, where 6 is 
the scattering angle of the electron with respect to the electron beam. 



that we can measure the angular distribution of the individual leptons. This might be 
useful because the leptons in Z-fusion tend to be forward, while the leptons from the 



signal are isotropic. This is illustrated in Fig. |T2|, where we show the angular distribution 
of the electron with respect to the electron beam for the same parameters as in Fig. 0. 
We see that 6 can be a decent discriminator between signal and background. 

And finally, given the current limits on rrih, the 6-fermion final state is not observable 
at LEP. From Table ^ we estimate that we can record 5 signal events at the NLC with 
^/s = 500 GeV and total integrated luminosity 50 fb~\ only if At < 300 GeV. Such low 
values of A^- are already ruled out by LEP (see eq. ( |4.24| ) above). 



4.2 Hadron colliders 

At hadron colliders we can similarly have 2-jet, 4-jet and 6-jet final states, by replacing 
e — i> g in Figs. ^ |] and |^. 

The dijet signal should be studied analogously to the graviton exchange |^, the 



difference now being that the angular distributions should be characteristic of spin 
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rather than spin 2 virtual exchange. 

The observation of the 4-jet final state will be rather difficult. There will be large 
irreducible backgrounds from tt, Z + 2j, WZ and ZZ, besides, the signal looks just like 
gg ^ h with extra radiated jets. In any case, the presence of extra jets is not very useful 
at a hadron machine. Even for larger Higgs masses, where the Higgs decays to gauge 
bosons, our signal is not distinctive enough by itself and will be simply added to the 
samples for the regular Higgs search. 

The 6-jet final state is cleaner, but the cross-section is too small to be relevant for 
the Tevatron. However, Higgs pair-production processes have small Standard Model back- 
grounds. Therefore identification of both Higgs particles in the event can be an important 
test for the presence of some new physics (see, e.g. 1^^). The main Higgs pair-production 
process at the LHC is gluon fusion [Q. In the four b-jet channel, very high pt cuts on the 
jets are necessary in order to suppress the QCD backgrounds, and then the acceptance 
is only 1-2%, assuming rather optimistic b-tagging efficiencies [0]. Therefore, for low 
Higgs masses it seems rather difficult to see a signal, and almost impossible to distinguish 
it from hbb at large tan (3. For heavy Higgs masses, where the cross-section is smaller 
than the one shown in Table we have 4 gauge bosons + 2 jet signatures. Then, for 
the purposes of both triggering and reducing the background, one would have to use the 
leptonic decays of a sufficient number of W^s and Z's, thus paying a significant price in 
branching ratios. 



5 Conclusions 



The possibility that new large compact dimensions may be probed in the near-term future 
at collider experiments is very exciting. Finding KK excitations will provide a strong 
evidence for the existence of extra dimensions. On general grounds, one expects that in 
addition to gravitons, there are other particles feeling the extra dimensions. If they are 
some SM fields, then looking for their KK excitations will provide the strongest probe 
of the extra dimensions. The experimental constraints on the masses of the KK modes 
of the SM fields require the corresponding extra dimensions to have sizes smaller than 
~ TeV~^. If we assume that all SM gauge fields are localized on our brane, any bulk 
field must be neutral under the SM gauge group. Hence bulk fields will only couple to 
some neutral combination of the SM fields. In the case of a bulk fermion, it may couple 
to iiH and be identified as right-handed neutrino. The smallness of the neutrino masses 
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can then be explained by the large size of the extra dimensions [BR B3] . It may also give 



rise to invisible Higgs decays [|l^]. If the bulk field is a scalar, it can couple to neutral 
combinations of the SM fields such as fifRH, F^^ F^^, or higher dimensional ones QQQL 
and so on. The study of this paper can apply to any scalar bulk fields which couple to 
JlIrH, independent of whether there is a flavor theory. The scalars which couple to 
F^'^F^y are similar to the dilatons and will be left for future studies. The scalars which 
couple to higher dimensional SM operators like QQQL, have suppressed interactions and 
are therefore much more difficult to produce in collider experiments. 

The new signals studied in this paper are not limited to the ideas of large extra di- 
mensions. Ignoring for the moment the difficulties in flavor model building, the results 
are also applicable in the usual four dimensional flavor theories if the flavons are light. 
In spite of the differences in the invariant mass distribution of the flavon (or its decay 
products), it seems that distinguishing between a four dimensional and a higher dimen- 
sional flavon scenario is intrinsically quite difficult and may have to wait for even more 
advanced experimental facilities. 

To summarize, we have studied many new signals coming from flavor physics within 
the framework of large extra dimensions and TeV scale gravity. These new signatures are 
particularly interesting because they are associated with the Higgs boson, which is the 
only missing piece of the Standard Model and hence the most hunted particle. We find 
that the production of Higgs + real flavon could probe the string scale just as well as the 
KK graviton production at the current and future colliders, if the flavons and the graviton 
live in the same number of extra dimesions. However, if the flavons only live in a subspace 
of the extra dimensions accesible to the graviton, associated flavon-Higgs production can 
provide stronger constraints than graviton production. In addition, there is a range of the 
parameter space in which the flavons can have macroscopic decay lengths, which will give 
rise to the spectacular signals of displaced vertices. The virtual flavon exchange generates 
higher dimensional operators which contain Higgs flelds. However, they do not provide 
as good probes as the real flavon emission, because the backgrounds are very large for 
flnal states with or 1 Higgs bosons, while the clean signals of 2 Higgs flnal states are 
too small to be observed. A possible exception is the case in which the flavons only live 
in one extra dimension. The fermion pair invariant mass peaks at the mass of the flavon 
zero mode, providing a very distinctive signal. This is a novel feature of the flavons in 
one extra dimension, since one extra dimension for the graviton in the TeV scale gravity 
scenario is ruled out. Note that in our study we have made the conservative assumption 
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that each flavon only couples to one corresponding fermion pair. There could be more 
promising signals if one relaxes this assumption. For instance, if a flavon couples to both 
quarks and leptons, then at a hadron collider we can have bb (from h) + leptons final 
states, which is an exceptionally clean signal. 
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A Formalism and conventions 



Expanding H and x around their vevs, H = (0, {v + h)/^/2) (setting the would-be- 
Goldstone fields to zero) and Xn = (Xn) + Xn; ^^e interactions ( |1.7| ) contain 

Hint D mfhfR + ^/^/^^ + E -j£fj^f^^^ + E Ti^^^^^^^n + h-^-' (A-25) 

where f = 246 GeV is the Higgs vev, mj and \f are the mass and Yukawa couphng for 
the corresponding fermion, h is the real Higgs field and x!n complex KK states. For 
simplicity, we assume that the masses of the real and the imaginary parts of x' are the 
same. The mass for the KK state Xn is given by = rn^ + 4'7r^|np/L^, where mg is 
the mass of the zero mode. Also, for notational simplicity, we will omit the prime on x' 
for the most part of the paper: x is understood as the vev shifted field when there is no 
confusion. 



We follow the approach of Han, Lykken, and Zhang changing slightly their no- 
tation: R ^ L for the length of the extra dimensions and n — 5 for the number of extra 
dimensions accessible to x's- Because there is a mass uiq for the zero mode, the number 
of KK states of x i^i the mass interval rfm^ is now given by 

dN ^ p{mn)dml, (A.26) 

with 

rS(^2 _ 2\(S-2)/2 

The effective interaction due to the virtual x KK state exchange can be obtained by 
summing over the KK state propagator, in a way similar to the virtual KK graviton 
exchange |^, ^ For 6 > 2 the virtual exchange is dominated by the high KK states, 
it depends sensitively on the cutoff scale Ac, but not on mo, as long as mo is much smaller 
than Ac (and ^/s, ^J\t\ for 6 = 2). The sum over the KK mode propagators gives 

Dis) = dml p{m^) ^— (A.28) 



TT — z In ( " I 



s—mT, 



2L'^(s-mg)('^-2)/2 
(5-2)r(<5/2)(47r)^/2 



IT -t\ ^ 

> s—mf' 



5 = 2; 

A._^A(^-2)/2- 



6>2. 



The real part represents the real KK mode production. It contributes to the same pro- 
cesses as the virtual exchange only if the KK mode of x has prompt decay. It is smaller 
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than the virtual contribution if ^ ^/s, rriQ for 5 > 2. The t-channel exchange can be 
obtained similarly by replacing s by t (remembering that t < 0). 

The virtual exchange for 5 > 2 is independent of s and therefore can be written as an 
effective dimension 8 operator, assuming Ac ^ rriQ and using Mp = M'^^^L^ , 

{5 - 2)r((5/2)(47r)V2 m^+^ '^^^"^ ' 
= -jjjhfRhW (A.29) 

Experiments on virtual exchange only probe the coefficient of this effective operator, A^, 
which is a combination of the gravitation scale M and the cutoff A^. Because the cutoff is 
not known, the experimental probes of virtual exchange processes should not be directly 
compared to the real x production. 

For 5 = 2, the virtual x exchange still depends logarithmically on the momentum flow 
through X- In analogy to ( |A.29| ), we will write the effective interaction as 



\yf\ 



2 



In 



— rriQ 



\lLfRh\\ (A.30) 



SvrM^ 

where is the momentum flow of the virtual x- 

For 5 = 1, the integral over KK propagators is convergent, 

D{s) = ^^^= + 1-^^ ■ (A.31) 

^ - tuq ^ y A^ - mg 

If Ac is taken to be inflnity, the second term can be dropped. For s > mg, it is dominated 
by real x production. The flnal state depends on whether x decay inside or ouside the 
detector. For s < or t-channel exchange, the effective interaction is 

\yf? 1 



4M3./^2_p2 



\fLfRh\\ (A.32) 



The singularity at = ttiq will be smoothed out by the decay width of the x KK states. 

There are also interactions with one or both h replaced by their vev. For 6 > 2, the 
effective operators ( |A.29|) are independent of the momentum. They can be implemented 
in COMPHEP via the exchange of an auxiliary fleld, without the need to consider t- 
channel and w-channel contributions separately. They are automatically accounted for 
by COMPHEP in the process of constructing all Feynman diagrams for the particular 
process. 
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